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In this paper we investigate the vacuum energies of several models of quantum fields interacting 
with static external currents (linear couplings) concentrated along parallel branes with an arbitrary 
number of codimensions. We show that we can simulate the presence of static charges distributions 
as well as the presence of classical static dipoles in any dimension for massive and massless fields. 
We also show that we can produce confining potentials with massless self interacting scalar fields as 
well as long range anisotropic potentials. 



I. INTRODUCTION 



It is well known in the literature that the coupling of a bosonic field (spin or spin 1) with two static external 
currents with the shape of Dirac's delta functions concentrated at specific points of space gives rise to the Coulomb 
potential for massless fields, and the Yukawa potential for massive fields ll,y,y]- This fact is one of the great triumphs 
of the Quantum Field Theory and enables us to comprehend, at the quantum level, the interaction between classical 
objects by describing them with Dirac's delta functions. 

Some natural questions that arise from these results concern the coupling between quantum fields and external static 
currents concentrated at specific regions of space. We can mention, for instance, what might be the physical meaning 
of a spinorial static current with the shape of a Dirac's delta function, or which kind of systems can be described by 
the interaction between bosonic fields with other types of external currents spatially concentrated, instead of currents 
with the shape of a Dirac's delta function. Even for the simple, and well known, system composed by two point-like 
Dirac's delta functions in 3 + 1 dimensions, which leads to the Yukawa potential, we can ask about the influence due 
to additional typical quantum interactions, like the A^** model in the case of scalar field. 

In this paper, which is devoted to the study of some of these questions, we calculate the vacuum energies for several 
models of quantum fields in interaction with external static currents concentrated at distinct parallel branes with 
' arbitrary number of codimensions. Specifically, we have treated the cases in which static delta like currents interact 
linearly with massive or massless fields. We would like to point out that systems of quantum fields interacting with 
external potentials concentrated along branes (quadratic couplings in the quantum fields) have been treated in the 
literature from time to time and in different contexts, see for instance P, 0, and references cited therein, but, 
as far as the authors know, the coupling of quantum fields with external currents concentrated along branes (linear 
couplings in the quantum fields) is not a well explored subject. 

The paper is structured as follows; in section (|Tl]) we calculate the vacuum energy of a real scalar field with mass in 
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d + D + 1 dimensions and coupled to N external Dirac-like currents concentrated at distinct regions of space. That 
is a generalization of some results which can be found in the literature for 3 + 1 dimensions 1] and are important 
calculations also because we establish, in a simple system, notation and the mathematical tools used in the whole 
paper. In section (jllip we consider the scalar field interacting with an external current with the shape of Dirac's delta 
functions derivatives. The results obtained are interesting once they show that the model considered can be used to 
describe the interaction between electric static dipoles in d + D + 1 dimensions. Section pV[) is devoted to extend 
the previous results to the electromagnetic field, that brings the previous results to more realistic contexts. In section 
(|V)) we consider a model in d + D + 1 dimensions described by the fermionic field interacting with a static external 
current with the shape of Dirac's delta hmctions. We show that, for the massless field in 3 + 1 dimensions and a 
current composed by two point-like Dirac's delta functions, the vacuum energy exhibits a similar spatial behavior in 
comparing with the electric field produced by a static electric dipole. In section (jVip we investigate the corrections 
due to the Xcf)'^ self-interaction on the results obtained in sections (|TT]) and piip for the restricted case with 3 + 1 
dimensions and currents concentrated along two branes. We show that we can produce confining potentials in first 
order in A and, as a particular case, the correction to the coulombian potential is proportional to the distance between 
the charges. Section (jVII[) is devoted for some final comments and conclusions. 

Along the paper we shall deal with models in d + D + 1 dimensions and use Minkowski coordinates with the diagonal 
metric (1, —1, —1, — 1). 

The time coordinate shall be represented by and the four-vector position shall be designated by 

a;= (a;°,x\...,x^x'*+\...,x'^+^) . (1) 

We shall also use the following notations 

= , (2) 

and similar ones for the momentums k. 

II. SCALAR FIELD IN d + D + 1 DIMENSIONS WITH DELTA-LIKE CURRENTS: CHARGES 

DISTRIBUTIONS 

In this section we study the vacuum energy of a system composed by the real scalar field in D + d + 1 dimensions 
with mass and coupled linearly with N d-dimensional delta functions concentrated at different regions of space. The 
model is described by the Lagrangian density 

where we have N fixed d-dimensional spatial vectors a^, p = 1, and xj_ is defined in ([2]). The parameters dp are 
the coupling constants between the field and the delta functions. 



Notice that the model ([3]) can be interpreted as the scalar field interacting with the external current 

N 

Ji{x) = ^cTp (5'^(x_L - Sip) , (4) 
p=i 

(or a sum of external currents) so the functional generator of the Green's functions can be written as 

Z, = exp(-^ / / d-^+^+^x d^+^+iy Jiix)Ap{x,y)Ji{y)\ , (5) 



where Ap{x,y) is the Feynman propagator 



ik{x — y) 



Af(x, y) = I (2^^pi5TT p _ ^2 ■ (6) 



^d+D+iy. exp 

.2 _ ^2 

The functional generator of any quantum system, at the limit T — )■ 00, can be written in the form 0] 

Z = eyiY>{-iET) , (7) 

where E is the lowest energy of the system. 

Comparing Eq's ([5|) and ([7]), this last one with E — Ej and Z — Zj, we have 

^' = tZ^ J J d'+''^'^d'^''^'y Ji{x)Ap{x.y)Ji{y) . (8) 
Substituting (O and (g]) into ^ leads to 

N N ^ 

El ^^^(yp<yq ^^P,q , (9) 
p=l q=l 

where Xp^q is the following integral 

Xp,, = J J d°+'+^x d^+'+^y 5^"^ (x^ - ap) (x, y)S<^''^ (y^ - a,) . (10) 

Expression ([9]) contains terms where p — q which account for contributions for the energy due to the self interactions 
of each delta function with itself. These terms must be discarded once, each of them, can be interpreted as the self 
energy of a given brane and do not contribute to the force acting on any brane [10]. So the energy due, strictly, to 
the interaction between the deltas is given by 

N N ^ N ^ 

p—1 q—1 p—1 

N N ^ 

= ^^(^pf^qO^- ^pq)l^'lp:q ^ (H) 
p=l q=l 

where 5pq is the Kronecker delta. 



In the appendix A the integral defined in pop is simplified to the form 



d^k 1 

Ip.g = -T / j^^^^2 exp (-jki • ap,) , (12) 



where we defined 

Substituting the result (fT2|) in Eq. (|lip we can write 

^/ = § = -^I]I]^P^g(l-'5p,) expHk^-ap,) , (14) 

p— 1 g— 1 ^ 

where we denoted by £/ the vacuum energy per unit parallel volume — J d^yn ■ 

From now on it shall be convenient to study, separately, the cases with and without mass. For the massive case the 

integral which appears in Eq. (jl4p is ill defined for higher dimensions, so wc shall search for its analytic extension in 

such a way to obtain a well defined expression for the energy (jl4p valid for higher dimensions as well. For this task 

we use the results of appendix B ([55)1 and ([M]) . this last one with ly = {d/2) — 1 and p = k^, in order to write 

= {2T:f/^m'^'^{may-'-'^/^^K(^d/2)-i{ma) , 0<d<5, m>0, 

= {2TTf/^m'^-^Gd{ma) , 0<d<5, m>0, (15) 

where Kn{x) designates the ii'-Bessel function and we defined 

Gd(a;)=xi-W2) i^(,/2)_i(x) . (16) 

Notice that, even though the result (jlSp had been obtained with the restriction < d < 5, the right hand side of 
the last line of (jlSp is well defined for any integer d, so it is an analytic extension, valid for any integer d, for the 
integral on the left hand side of the first line[ll| in such a way that 

/ d'^ki — ^exp(±ik_L • a) ^ {2TTfl'^ m'^~'^Gd{ma) , to > . (17) 
J ^1+ 

Defining the variables apq — \sbpq\ = ja^ — a^j and substituting the integral which appears in p4p by its analytic 
extension pT|) the energy pi)) reads 

£i{ap, ag, Qpq, m y^O, d,N) ^ ---——^^^apaq(l- 5pq)Gd{rnapq) , (18) 

' p=l 9=1 

which is valid for any integer d. 

The force per unit parallel volume acting on a given current localized, for instance, by the perpendicular vector 
af can be calculated from expression (fT8|) as follows 

( d ^ \ „ 



Ca^fc ^ / agk 

N 



■ E ^ (™)^/^if./2(™) — , (19) 



where we used the fact that 

= (20) 

and discarded a term (1 — Ski) once the summation does not include the terms where k — i. It is interesting to notice 
that the force satisfies the superposition principle for the delta-like currents. 

The energy for the massless case could be obtained by taking m = in Eq (fTS]) . which is well defined for d ^ 2, as 
shall be exposed later. For complentess we shall perform the calculations with m ~ following a different approach. 

For the massless case let us consider, separately, the situations where d = 2 and d ^ 2. When d 2 we use the 
results obtained in appendix B (IM)) and (jlOOp . this last one with /i ~ id/2) — 2, v — {d/2) — 1 and p = k^, in order 
to write the integral which appears in Eq. (jl4p . for m = 0, in the form 

y" d'^ki exp(±zki • a) = i27Tf/^a^-'^ duu'^'^^^^-^J^d/2)-i{u) 

= (27r)'^/2 2(<i/2)-2 r - 1^ a^-d ^ 2 < d < 5 , (21) 

Similarly to what happened for (fT5|) . the result (pij) was obtained with the restriction 2 < d < 5, but the result in 
the second line of (|2ip is valid for d 7^ 2,0, —2, —4, so it is the analytic extension for the integral on the left hand 
side of (I^TI) and we can write 

1 / ^ \ 

(22) 



J j^exp(±iki • a) [2T,fl^ 2W2)-2 r|^^ - 1^ c?'" , d 7^ 2 , 

where we considered only d > 1. 

Substituting the integral which appears in (|14p with m = by its analytic extension (I22|) we have the energy 



r((d/2)- i)2W2)-3 

f7(ap,a„ap„m = 0,d^2,iV) = ^ ^--^^ E E {^-^v<i) ■ (23) 



)-J AT N 



0=1 q=l "-PI 



Before we study the case where m — and d = 2, let us investigate the behavior of the analytic extension (|17p for 



any d > 1 and m = 0. For this task we have to use the fact that 



i 



.=0 r(^)2^-i 

K^{z) — > , ly^O 

Koiz) ^ -ln(^|)-7 , (24) 

where 7 is the Euler constant. 

Substituting the first Eq. in ITBl) and pT)) we are taken to (1^^ . Similarly, when m — the energy is 

reduced to Eq. ([^5]) . So the result ^5)1 can be considered valid also when m = and d 7^ 2. 

For the situation where m = and d — 2 the energy (|18p exhibits divergences which do not depend on the distance 
between the branes and can be removed with no relevance for dynamics of the system. For this task we substitute 



the second Eq. (I24p in and (|T8l) what leads to 
£i{ap,aq,apq,m-f 0,d = 2,N) = ■^^^<yp(7g{l-Spq)liii j +^^^crpaq{l-5pq) 



InR^ +7 



, (25) 



where oq is a non vanishing arbitrary constant with dimension of length. 

In the limit to — > the last summation on the right hand side of (j25p is composed by divergent Cpa-independent 
terms, which do not contribute to interacting forces, so it can be discarded whith no relevance for dynamics, 



N N / \ 

p=l g=l \ / 



£i{crp,(Jq,apg,m = 0,d= 2,N) ^ —2_^2^apaq{l - Spg)lii[ ^ ] . (26) 

p=l q= 

The calculations for the massless case with d — 2 can be interpreted as an insertion of a mass as a regulator 
parameter introduced in order to remove a-independent divergences. 

The force in the massless cases is obtained following the same steps which lead to the result (fTOj) . So, from 
and (|^ we obtain 



r(d/2)2W' 



2)~1 



N 



which is valid for any value of d. 

As a check of consistence we can take the limit of vanishing mass in Eq. by using Eq's ([M)) . The results are 
given by the force ([27|l for any value of d. 

The situations with 3+1 dimensions and N = 2 are of special interest; when d = 3, D — and iV = 2 we have two 
point-like currents in 3 + 1 dimensions, whose interaction energy is the Yukawa potential which can be obtained from 

m B 

Ej{a, m,d = 3,D = 0,N = 2)^ _ai^ exp(-TOa) 

in a 

where we suppressed the sub- indices 1, 2 for a, once there is only one distance involved. For m = the above result 
leads to the Coulombian interaction. 

For d = 2, _D = 1 and iV = 2 we have two delta-like currents concentrated along two different parallel strings placed 
at a distance a from each other. In this case the energy per string length reads 

£i{a,m,d^2,D ^ 1,N = 2) = -^-^Ko{ma) , (29) 

27r 

which is reduced, in the case m = 0, to the expression 

£:/(a,TO = 0,d = 2,i:)=l,iV = 2) = ^i^ln( — | , (30) 

2tt \aoJ 

where we discarded a divergent a-independent and oq is a non vanishing arbitrary constant, which does not contribute 
to the force. 
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The case where d = 1 , D = 2 and N = 2 corresponds to the situation of two delta currents concentrated on parallel 
planes. In this situation the energy per unit of area (jlSp becomes 



£j{a,m,d^l,D ^2,N = 2) = -^^exp(-ma), (31) 

2m 



which reads for vanishing mass 



8i{a,m = 0,d=l,D = 2) = , (32) 



after discarding a divergent a-independent term. 

III. SCALAR FIELD IN d + D + 1 DIMENSIONS WITH DELTA-LIKE DERIVATIVE CURRENTS: 

DIPOLE DISTRIBUTIONS 

In this section we study the model described by the Lagrangian density 

Ch = - \m'<t>^ + {^''Y/'v ^^Krt '^'(^i - ^v)^'^ - (33) 

where V^^^ designates a fixed and static four vector in the reference frame we are performing the calculations. 
The Lagrangian p3|) can be interpreted as the scalar field coupled to the current 

N 

Jii{x) =d^J2''p ^Mp) '^''(^-1- - ap) ■ (34) 
p=i 

Following similar arguments which led to equation (jlip and performing two integrations by parts we have that the 
energy of the system (|33)) is given by 

N N ^ 

= X! X! '^P^'? (1 - ^pq)7p^J'pq ' (35) 



2T^ 

P=l 9=1 



where jTpg is the integral defined by 



Jp, - / / d^'+'+^x d^'+'+'y <5('^)(xx - s.p)(vi;^Vl'^)d,d. AF{x,y)y^^Hy± - a,) • (36) 



In the appendix A it is shown that 

Jpq = {^{p)± ■ ^pqJ-) (^(g)-L • "^pqA^lpq , (37) 



where V(,)^ - (V^^^^, l^^^^^) and 



^P9i - ( TJTT'-'lJ^r ) ' (38) 



with a' designating the z-th component of the vector a.p 



By substituting Eq. dST]) in and with the aid of ([H]), and HH) we have finally 



N N 



^11 = ^££crpCr9(l-'5p<?)(V(p)± • Vp,;i)(V(q)^ • V 



P=l 9=1 



m G{mapq) 



1 



■^^crpcrq(l - 5p 



(27r)'i/2 



^ ^ p—1 q—1 

-(map5)"^"^'^/^'ifi+(d/2)(TOapq)(^V(p)_L • (map,)^ ' ( 



■P9, 



(39) 



For the massless case and d ^ 2 we proceed in a similar way as we have done to obtain the energy p9p but, instead 
of pT|). we use the result ((^ what leads to 



2(d/2 



W AT 



-r(d/2)^^apa,(l-5p,)-^ 



p=l q=l 



V 



{p)± ■ "pij 



V 



(q)± ■ a-pg 



d 7^ 2 . (40) 



The situation where m — and rf = 2 is obtained by taking the limit m — > of (j39p with d = 2, what is done with 
the aid of ([M)) . as follows 



1 ^ ^ / 

" '^P?) (V(P)-L • ^P9-l) (V(9)± • Vpq_L) In ^ 
p=lq=l V 



p=i q 



P9 

ao 



47r 



^^o-pCrg(l-(5pq)- 



p=l g=l 



pq 



V 



(p)± • Spg 



V 



■ ap5 



d = 2 , 



(41) 



where, similarly to what was done in the previous section, in the first line of the above equation we discarded an 
flpg-independent divergent term and introduced an arbitrary finite constant aq with length dimension. 

Notice that if we take d = 2 in Eq. P0|) we obtain it is, Eq. PD|) can be considered to be valid with no 

restrictions for d. 

In order to have a better insight on the meaning of the current (IM)) . as well as of the energies (1591) and ([3D]), let us 
define a — ai2, designate by a the unitary vector in the direction of a and consider the restricted situation for Eq. 
(IM)) where = 2 and d = 3, 

exp(— ma) 



Siiid = 3,7V = 2) 



(ma)2 + 3(ma + 1) ( ai V^^' • a ) ( a2V'f' ■ d] - {ma + l)[ aiV^^' ■ da V 



(2) 



r(2) 



An interesting situation occurs for the massless case, where equation (|42p reads 

1 



(o'iV(i) • a2V(2)^ - 3(^(TiV(i) • d^ (^(T2V 



(2) •« 



.(42) 



(43) 



Notice that Eq. (|43|) is exactly the interaction energy, with the negative sign, between two electric dipoles whose 
electric dipole moments are given by criV(]^) and (J2^{2)- 

The energy ([39]) is the generalization of the interaction between the dipoles distributions of the scalar field along 
N _D-dimensional parallel branes. 

It is worth mentioning that one could calculate the interaction between a dipole and a point-like charge by simulating 
them with currents of the form ([M]) and ([5]), respectively. 
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IV. MAXWELL FIELD 



In this section we extend the models exposed in sections (|ll]) and (jllip to the MaxweU field, trying to be as general 
as possible. 

The first model we consider is described by the Maxwell field Lagrangian interacting with an external current, as 
follows 



N 



(44) 



where is the photon field, F^^ — — and W^^p) is a four vector satisfying the condition = which 



assures that 



cTp Wy^(p) (5'^(xx — Sp) = and so, the gauge invariance of the Lagrangian (P^ . The quantities W^j 
shall be taken as constants and uniforms four vectors in the reference frame where the calculations are performed. 

Following similar steps of sections ([n| and (jllip , and using the facts that = and the photon propagator is 
given by 



fc2 



^'^'^-(1-7)- 



A:2 



exp[-ifc(a; - y)\ , 



we can show that 



p 2{d/2)-l f fj \ ^ { \ / \ 



(45) 



(46) 



For the case where d = 2 we insert a mass parameter in the propagator (|45[) as a regulator parameter, as follows 



It leads to the energy 



^r-(i-7)^ 



exp[— i/c(x — y)] 



(47) 



£ttt — 



E 



N N 



III 



(48) 



where, in the second line, we used the second equation (|24p . discarded an apa-independent divergent term and ao is 
a finite arbitrary constant. 

By setting N = 2, d^i, D = Q and = rj^'^ in Eq. (|46)) we have the Coulombian interaction between two 
charges in 3 + 1 dimensions, whose energy is 

(71(72 1 



E 



Att a 



(49) 



with a designating the distance between the charges. 
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The second model we shall consider for the Maxwell field is an extension of the one considered in section (IIII|) for 
the scalar field. Its Lagrangian density is given by 

11^ r n 

Hiv = —.F^^.F^^" - ^(a^^^)' + 5] ap iy^(rt a„ \8\^^ - ap)J . (50) 

Similarly to what happens for the model (l44l) . the same condition = assures the gauge invariance of (|50p . 
The quantities V^^^ are four vector with no restrictions and, as well as W^t(p) , are taken to be static and uniform in 
the reference frame where the calculations are performed. 

By similar procedures which lead to results (|¥T|) . (1^5]) and (|^5)) . the Lagrangian (|50p leads to the energy per 

unit of parallel area 

2W2)-2 ^ ^ 



p=l q=l 



1 



PI 



(51) 



which is valid for any value of d. 

The result ([5T|) can be interpreted as the interaction energy between static electric dipoles distributions along N 
Z?-dimensional parallel branes. This fact can be seen by considering the most interesting situation, which occurs for 

= 0, d 3, AT = 2 and T^^'^^ = M^^'^^ = 77^°, where equation ^ reads 



E 



IV 



1 1 

47r 



(52) 



which is the interaction energy between two electric dipoles in 3 + 1 dimensions separated by the vector a and with 
electric dipole moments given by cfi'V^i) and (72V(2). 

By these means each term in the summation present in (j50p can be used to simulate the presence of a static electric 
dipole distribution along a brane with D dimensions. 

V. DIRAC FIELD 



In this section we consider the interaction of the Dirac's field, ip, in d + D + I dimensions with an external static 
current concentrated at some regions of the space. The Lagrangian density of the model is given by 



Cv = Hh^d^ ~ m)ij + (7p ^ip)6'^{x± - ap) j + ^P cr, ^iq)S'^i^± - a,) j , 



(53) 



\p=l / \q 

where ^(p) are N fixed and static spinors in our referential frame, 7^ are the Dirac matrices and, as usual, ^(p) = ^(p)^ 
Notice that the Lagrangian (|53p can be interpreted as the Dirac's field coupled to the external fermionic current 

N 



X^'^P C(p)'5''(xi - ap 



(54) 



P=i 
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Using the fact that the functional generator of the model ([3]) is given by 
Z = exp 



rfd+D+i^ /^^^ -ap) e(,)'5'(x± - a,) 

/ \9=1 



(55) 



where 



Sfix,y) - (^7^a^ + m) / ^^-p^-^^^— — ^ exp Hfc(a; - y)] 



■ exp [—ik[x — yj\ 



(56) 



is the usual Dirac propagator, and proceeding as in the previous cases, it can be shown that the field energy per unit 
of parallel area is given by 

N N 



= 1^ = ^ (1 - ,5p,) ^(p) / 7^ ^2^^+ t*'^^ • M 

p—1 q—1 L-/ V / ± 



Defining the spatial vector 



*Up, = e(p)7 7 = (7\---,7''+^) 



(57) 



(58) 



equation (|57|) becomes 



p=l q=l 



Upq • VpqJ_ - ™C(p)C(g) 



1 



exp [ik_L • (apg)] 



N N 



Up, • Vp,_L - 'mi(p)^(q) G{mapq) , 



p=l q=l 

where, in the second line, we used the result (|17p and definition ([1 
With the aid of (I16p equation ([5^ reads 

= -(27r)-'^/W-i^^(l-5p,)(map,)i-W^) 

p=l 9=1 

+-f(d/2)-l(™apq) (cr*f(p))(crg^(g)) 



i^£j/2(map5)l (CTpUpgCTg) • ^ 



Taking into account that ap, = -a^p , Op, = a,p, (a-p^(p)) (0-9^(5)) = {<^*qi(q)){(^pi(p)) and (cr*UpgCrg) 
'(cr*Ugp(7p)*, equation ()60|) can be written in the form 



(59) 



(60) 



N N 



^ ^ p—1 q—1 



Nl-(<i/2) 



Kd/2{mapq) 3 (cr;C(p))7 (o-9^(«)) 



+-K'(d/2)-i(»7iapg) 3? (ct*C(p)) (0-9^(9)) 



(61) 



where we used definition (|58p . 
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For the massless case with d = 3, D = and N = 2 the energy ([61]) reads 

1 1 a 



£5(771 = 0, d = 3) = £5(777 = 0, d = 3) 



27r a 



KC(1))7 (0-2^(2 



where a = 012, which leads to the force acting on the current 2 



2tt a3 



3d 



a ■ ^J(KC(i))7 (^2^(2))) - 3((o-t?(i))7 (^2^(2)) 



(62) 



(63) 



It is interesting to notice that the force field has the same spacial behavior as the one exhibited by the electric 
field produced by a static electric dipole given by 3^((t5'^(i))7 ((72^(2))^ • 



VI. Afli* MODEL WITH TWO DELTA-LIKE CURRENTS: CONFINING POTENTIALS 



In this section we consider the radiative corrections induced by the, well known, self-interaction on the energy 
of the scalar field interacting with static external currents in 3 + 1 dimensions. 
The Lagrangian of the model is written, in the general form. 



(64) 



where ^m^ is a mass counter term and the current depends only on the perpendicular coordinates x_l. Once we are 
restricted to 3 + 1 dimensions, we must always have d + D = 3. 

By using standard perturbative methods of quantum field theory we write the functional generator Z(3) up to first 
order in the coupling constant A 



Z(J) — exp 



i d X 



A 5" 



*5^P ~o / d X d y J(x)AF(x,y)J(y) 



(65) 



2 W(x)2 ; V4! 5J(x)^ 

From now on we shall consider the restricted case where the current J is located at two distinct regions of the space 
Substituting the functional (j65p in Eq. ([7]) we have, at first order in A, 



E ^ lim — / d' 

T^oo 2T 



d+D+l^ fld+D+l 



y J(x)Ap(x - y)J(y) 



^'^+'^+^x d'^+^+^y d'^+^+^z J(x)Af(x - y)AF(y - z)J(z) 



(66) 



where we disregarded non quadratic terms in J in order to avoid contributions to the energy due to self interactions, 
as shall be exposed later. 

The first term on the right hand side of Eq. (|66p gives the contribution to the energy in zero order in A studied in 
sections (jnl and (jlllj) . The second term gives the correction in order A to the energy and, by the use of ([6]) and after 
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the integrations d'^+^+^x d^+^yn , can be simplified to 

In above expression, the counter-term 5iv? shall be defined for each specific situation considered later. 

Now, let us take an external current with in the form (|4j), but with only two Dirac's delta functions present, 

J(xi) = a^5''{^i_ - ai) + fT2(5'^(xi - as) . (68) 



Substituting (j68p in (j67[) . performing the integration in d'^y±_ d'^z± and discarding the self interacting terms between 
a given delta with itself we have 

.(1) E(^) fSm^ 6iX \ f d-^k^ exp(-ik^ ■ a) 

1 / (5™^ gj;^ \ ^^^2 9 /■ d'^k± exp(-ik^ • a) 



2 V 2 4! ^ ' J m dm J {2ttY k^ + 

where we defined a = ai — as . 

Contributions for the energy (|66p quadratic in J give rise only to self interacting terms for the deltas, and must be 
discarded. Using the result (fTT)) Eq. (|69|) gives, for the case with mass, 

-Kd/2{ina) + {d- 2)(ma)"^X(rf/2)_i(ma) (70) 
The counter-term (Sm^ is determinated by imposing that the above expression is finite, what is equivalent of taking 

+ — Aj.(0) =AK(m,d) (71) 



2 V 2 4! 

where «;(m, d) is a given finite parameter which must be determinated experimentally, can depend on m and d and 
whose dimension is mass square. So, Eq. fTO]) reads 



£^^'){m,d,D,a) = -XK{m, d)-^^^a^-'^{ma)^'^^^^-^ x 

--ft^<i/2(™a) + (d- 2)(TOa)"^iir(d/2)-i(»7ia) . (72) 

For d ~ 3, d = 2 and d = 1 we have, respectively 

£:W(m,d^3,Z? = 0,a) ^ A^c(™, d = 3)^^^^^^ 

477 m 

fW(m,d = 2,Z? = l,a) = A«(m, d = 2)^0^^!^ 

27r ma 
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S^^'' (m, d ~ 1, D ^ 2, a) = \K(m,d — l)—^-—^a'^{l + ma)- — exp{—ma) (73) 

2 (ma)'^ 

The situation where 771 = can be obtained by taking the hmit of vanishing mass in above results, what can be 
easily done for d ^ 2. In this paper we avoid the case m = and d — 2, where infrared divergences are problematic. 
For d 7^ 2 wc take m = in the first line of Eq. ((69|) and use the fact that 



yd'^plexp(±7pa) = {2T:f'^a^-^ duu'^^'^^-^J(a/2)-i{u) 



(27r)''/V-'^2W2)-4pM_2] , d^2 , (74) 



which is obtained in appendix B. So we have 



- ( ^ + f A.(0)) (g^«-^2(^/^)-r - 2) , . ^ 2 . (75) 



For d — ?> and d — 1 the counter-term dm? is determinate by ((7T|) and we have 



f(^^(m = 0,d = 3,15 = 0, a) = -AK(m = 0,d=3) 



CTicr2 
( 

47r 



£:(^)(m = 0,d= l,i:> = 2,a) = \K{m = Q,d=l)^^a^ . (76) 

6 

We would like to point out that, depending on the signals of the parameters k and k! , the expressions given by (|76|) 
can be confining potentials. 

Another interesting current is the one considered in section (jIII[) which we shall take, in this section, to be composed 
only by two terms, that is, we shall take = 2 in ([33| . 



J(xi) = aiF(^)a^[(5(xi - ai)] + cj2Vl'^^d^,[5{^i_ - aa)] . (77) 

Substituting (|77p in (I67p we have, after some manipulations similar to the ones done in the previous sections, 
^n-i f drri^ 6iA , „ „ , Z""" d''kj_ exp(— ik^ • a) 

= - (— + ir^-(o) j ^^'^^(^(^)- ■ ^")(^(^)- ■ I (2^ iA^m^y ■ 

The above integral is calculated in and (|70p . 

The general result for (|78p is not much suggestive. The most interesting situation occurs in the massless case for 
d = 3, where we must take the counter term 5m^ in such a way that 

1 / 5ni? 6iX . , A , , , 

2 [~Y' + ir^^^^V " ' ^ ^ 

where k is a finite constant which must be determinated experimentally. So Eq. (I78p leads to 



47r a 



V(i)± • V(2)i - (V(i)i • a) (V(2)i • a) (80) 



which is an anisotropic long range interaction energy. In fact, its angular dependence is the same one exhibited by 
the interaction between two static dipoles. 
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VII. CONCLUSIONS AND FINAL REMARKS 

In this paper we have made an investigation on the role of the interaction of quantum fields with external static 
currents concentrated along different parallel branes. Specifically we generalized to d+D+1 dimensions the interaction 
between bosonic fields (scalar and electromagnetic) with external currents composed by N d-dimensional Dirac delta 
functions and delta functions derivatives. For the scalar case we considered the situations with and without mass. 
We showed that we can describe charges distributions as well as static dipole distributions by the use of currents 
concentrated along parallel branes with arbitrary codimensions. 

We also had made the same study for the fermionic field and, in the massless case, we showed that we have a 
dipole-like interaction. 

We showed that the first order radiative corrections induced by the A(/>^ self interaction, for the case of scalar field, 
can lead to confining potentials between point-like currents, as well as a long range anisotropic interaction. 

As a last comment we would like to say that four-pole interactions can also be described by the use of currents 
concentrated along specific regions of space ^. 
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Appendix A 



In order to write the integral (|10p in the form (|T2l) we, first, use definition ([6]), what leads to 

Ip^g = [ (dPxdPy I [ d^'x^^d^'y^^ [ /rfVdV ^^'''(x±-ap) ,5('^)(yi-a,) 



X 



/fr/(^/ (^exp[zfc°(x°-yO)]expHkr(x||-y,|)] 
exp [-zk_L ■ (xj. - y_L)] 



and integrating in the variables dfc" and (ik|| , we have 

d'^kj_ 1 
(2^ k]_+m 



(81) 



Performing, in the following order, the integrals d'^x±, d'^y^, d^x and d^x\\^ using the fact that 

f dx 

5{k) - / ^ exp(zA:x) , (82) 



Tp,q = - I dy" I d^yii / ' exp [-ik^ • (a^ - a,)] (83) 
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As usually done in quantum field theory, we designate the integral in the time variable by T = J d'lf' . We also 
designate the integral in the parallel volume by 



(84) 



With these considerations, Eq. (|83p is written in the form ([2 

Another useful integral is the one defined in ([55)1 which, following similar steps used to obtain p^ . can be rewritten 
into the form 



•Jpq — TL 



D 



(27r)° 



(V(pU-kx)(V(,u-k, 



exp(— ikj^ • ap 



where we defined V(,)^ = (V^^^^ , V^^^ , V^^^). 
By using the operator ([38]) Eq. ([85)) becomes 



/ \ / \ f c^^k 1 



With the aid of 1^ Eq. §^ leads to ([57|) . 



Appendix B 



In this appendix we calculate the d-dimensional integral 



Id{a) = J d'^p /(p) exp(±ipa) 



where f{p) is any function of the modulus p of the vector p which satisfies the property f{p) = f{—p). 
For this task let us consider, first, the case where d > 2. 

It is appropriate to use spherical coordinates in d dimensions, in such a way that 

p^ = p cos(0i) 

p^ = p sin(0i) cos(02) 

p'^ — p sin(6'i) sin(02) cos(6'3) 



(85) 



(86) 



(87) 



/ ^ p sin(6'i)sin(6'2)---sin(6'd-2)cos(6'd-i) 
p"^ ^ p sin(6'i)sin(6'2) • • •sin(6'd_2)sin(6'd_i) , 

where < < tt, < (?2 < tJ", ••■7 < 9d-2 < ti" and < 6d-i < 2tt. In this case we have 

d-l 



d'^P = {pf-^ dp \[ sin'*-('+i)(0,) de. 



(88) 



1=1 
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It is also convenient to use a coordinate system in such a way that a = (a, 0, • • • , 0). So, the integral ((87l) becomes 



Idia) 



dpip) fip){ / dOi sm''--'{ei) exp ±ipacos(6ii; 
\Jo L 

^ sm''-\e2)]lj^ dOs sin''-4(03)j ••• n rf^d- 



2 sin'^-3(^rf_2) 



2tt 



dOd- 



By using the results 



a 



(iff sin ''{9) exp ±ixcos{9) 



(89) 



/ dO sin™(ff) = 
Jo 

equation ([89]) reads, after some manipulations, 



i/2 r[(l/2)(m+l)] 
r[(l/2)(m + 2)] 



TO = 1,2,3,. 



1 



, d>2 



where we made the change of integration variable u = ap. 
For d — 2 dimensions, usual spherical coordinates leads to 



h = I dpp f{p) / de exp 
/o Jo 



±ipa cos(ffi) 



Noticing that 



d9 exp 



±ixcos(6'i) = 27rJo(a:;) 



and performing the change of integration variable u = pa we have 



2tt f°° f u\ 

h ^ ^ du uJo{u)f\ - 
a'' Jo \a J 



For d = I dimensions we use the fact that f{p) — f{—p) in order to write 

/•OO 

h ^ I dp f{p)exp{±ipa) 



du f[ — \ cos(u) , 



where we also made the change u = pa. 
Using that 



COS(W) = \l -y/uJ_i/2iu) , 



we have 



iu)f 



(90) 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 



(97) 
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Combining the results ((9T|) . ((94|) and (|97|) we are taken to the general result 



d'^P f{p) exp(±ipa) 

Applying formula (|98p for the case where /(p) = (p^ + ni^)'^ and using the fact that 



du 



/o + x-^ 

which is valid for a; > and —l<iy< 3/2, we are taken to the result (fTS]). 
Taking /(p) = in Eq. ((98)) and using the fact that [8| 

/■°° , ^ r(l/2 + :y/2 + /i/2) 

/ du u^.Mu) = 2^ ^-Vt 

io r(l/2 + i./2-/./2) 

which is valid for — i^ — 1 </i< 1/2, we obtain Eq. (PTjl . 

Another case of interest corresponds to /(p) = p""* which, with the aid of (|100p . leads to ([7^ . 



(98) 



(99) 



(100) 
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